Noise in Bose Josephson junctions: decoherence and phase relaxation 
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Squeezed states and macroscopic superpositions of coherent states have been predicted to be 
generated dynamically in Bose Josephson junctions. We solve exactly the quantum dynamics of such 
a junction in the presence of a classical noise coupled to the population-imbalance number operator 
(phase noise), accounting for e.g. the experimentally relevant fluctuations of the magnetic field. We 
calculate the correction to the decay of the visibility induced by the noise in the non-markovian 
regime. Furthermore, we predict that such a noise induces an anomalous rate of decoherence among 
the components of the macroscopic superpositions, which is independent on the total number of 
atoms, leading to potential interferometric applications. 

PACS numbers: 03.75.Gg,03.75.Mn 



I. INTRODUCTION 

It has been realized in the last decade that an ultracold 
Bose gas trapped in an optical potential offers the pos- 
sibility to manipulate coherently entangled many-body 
quantum states, with interesting applications in preci- 
sion measurements and in quantum information. For in- 
stance, spin squeezed states and macroscopic superposi- 
tions of coherent states are generated by the dynamics 
of a Bose Josephson Junction (BJJ) [iMij]. The useful- 
ness of squeezed states in improving phase sensitivity in 
interferometry has been demonstrated in a recent exper- 
iment Q • An even better sensitivity is predicted to arise 
by employing macroscopic superpositions The pres- 
ence of noise and of coupling with the environment causes 
decoherence and limits the experimental time for coher- 
ent manipulations. Decoherence may even prevent the 
production of certain entangled states, a fundamental is- 
sue in the quantum-to-classical transition [6]. It is thus 
important to study the robustness of these nonclassical 
states in the presence of noise. Several sources of de- 
coherence in BJJ like particle losses |7j, collisions with 
thermal atoms [H, 0|, interaction with the electromag- 
netic field llOL and random fluctuations of the trapping 
potential [11[ have been identified and analyzed theoret- 
ically. 

In this work we solve the quantum dynamics of a BJJ 
in the presence of a noise coupling linearly to the number- 
imbalance operator. This noise results from the fluctua- 
tions of the optical potential and of the magnetic field, 
which are, together with atom losses, the main sources 
of decoherence in the experiments of Ref. 0, [H, [Hj. 
Our solution is exact and allows in particular to capture 
the decay of the Ramsey fringes visibility at short times 
(non-markovian regime). Furthermore, it shows that the 
macroscopic superpositions of phase states generated by 
the unitary dynamics are rather robust with respect to 
the noise considered. According to the usual scenario 
for decoherence [H], [Hj, by increasing the intensity of 



the noise these superpositions should be transformed into 
statistical mixtures of the same phase states at a noise in- 
tensity proportional to a negative power of the number of 
atoms, which characterizes here the "distance" between 
the phase states. For the aforementioned noise, we find 
that this is not the case. The typical noise intensity at 
which the coherences between the distinct phase states 
are lost is independent of the atom number and equal 
to (or even, for many-component superpositions, larger 
than) the noise intensity at which phase relaxation oc- 
curs. Phase relaxation means that each phase state of 
the superposition converges to a mixture of Fock states 
and acquires a completely undefined phase. At interme- 
diate noise, the phase has spread significantly but some 
entanglement remains due to the non vanishing coher- 
ences among the phase states of the superposition. We 
quantify this entanglement by computing the quantum 
Fisher information [5| and estimate the gain in phase 
sensitivity with respect to separable states. 

The paper is organized as follows. After introducing 
in Sec. II the two-mode approximation for the BJJ, we 
review in Sec. Ill the quenched dynamics of the junction 
in the absence of noise, leading to the formation of non- 
classical states. Sec. IV presents our results on the effect 
of the noise on the density matrix of the atoms and on 
the visibility of the Ramsey fringes, while Sec.V analyzes 
the degradation of the coherence of macroscopic super- 
positions. Finally Sec. VI offers a summary and some 
concluding remarks. 

II. MODEL 

We describe the BJJ by a two-mode Hamiltonian [l6| 
= X Jl ~ XJz ~ 1KJ X , (1) 

where the angular momenta operators J X1 J y , and J z 
are related to the annihilation operator aj of an atom 
in the mode j = 1,2 by J x = (a\a2 + a^ai)/^, J y = 
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i{a\a2 — a2tti)/2, and J z 



(a\di — ^2^2) /2, the at time t > reads 



latter being the number imbalance operator. We assume 
a fixed total number of atoms TV, i.e. we do not account 
for atom losses. We take N to be even for simplicity, the 
odd case being qualitatively similar. The Hamiltonian 
(PQ) models both a single-component Bose gas trapped in a 
double- well potential [l2[ - external Josephson junction - 
and a binary mixture of atoms in distinct hyperfine states 
trapped in a single well 0, [ItJ - internal Josephson junc- 
tion. In the external BJJ the two modes i correspond to 
the lowest-energy spatial modes in each well. For the in- 
ternal BJJ, the two relevant modes are the two hyperfine 
states. The first term in (pQ) describes the repulsive atom- 
atom interactions; for the external BJJ, x is the half of 
the sum of the interaction energies Ui in the two modes, 
whereas for the internal BJJ x — (Ui + U 2 )/2 — Uyi also 
depends on the inter-species interaction U12. In both 
cases, A = AE + (N — l)(U 2 -Ui)/2 is related to the 
difference AE = E2 — E\ between the energies of the two 
modes. The last term in (pQ) corresponds to tunelling be- 
tween the two wells or, in the internal BJJ, to a resonant 
laser field coupling the two hyperfine states. Both x an d 
K are experimentally tunable parameters. 

It is convenient to characterize a state of the BJJ by 
its Husimi function Q(0, 0) = ^|-0)| 2 on the classical 
phase space (the Bloch sphere of radius N/2), where 
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is a SU(2) coherent state [18], a = tan(0/2) exp(— 
and \n) is the Fock state satisfying J z \n) = n\n). 
For a coherent state ([2j) Q is peaked around the vec- 
tor N(smO cos </>, sin sin </>, — cos#)/2, the components of 
which are the expectation values of J x , J y , and J z in this 
state. In particular for a phase state, i.e. a coherent state 
with = 7r/2, such a peak is located on the equator of 
the Bloch sphere. A Fock state has a ^-independent dis- 
tribution, with a peak in 6 at = arccos(— 2n/N). 



^ \t) = ±(ji°\t))=co S N -i( X t). (3) 

At small times, the BJJ is in a squeezed state [lo[. Later 
on, it returns to the initial state \a = 1) at the revival 
time T = 27r/x; at intermediate times t q = T/(2q) it is 
in a superposition \^°\t q )} = uq Ylk=o Ck\e~ l27rk / q ) of q 
phase states @,[3| with \u \ 2 = 1/q and c k = e i7rfc ( fc + A/ ")/ ( ? 
(we have taken q even). By ([2]), the matrix elements of 
the density matrix p {0) (t q ) = |^ (0) (*</)> W> (0) (*</)! in the 



Fock basis are the sum over all fc, k' = 0, ■ 



1 of 



1 1 



(«(^)K) = -^v 



N 



s, —2iit(kn—k'n') /q iirik' 



N 

-U — 

^ 2 

~k' 2 )/q 



(4) 



with p£l(t q ) = q- 1 c k c% t \e- i2 ' Kk ^){e- i2 ' Kk '/% Since the 
dynamics does not couple the two modes, J z is a con- 
stant of motion. Thus (n\p( \t)\n) is constant in time 
and equal to P a =i(n) = 2 _Ar ( n+ ^ 2 )- OI> der to ad- 
dress later on the decoherence and phase relaxation of 
the superpositions of phase states, we decompose p(°\t q ) 
as 
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k=0 k^k'=0 
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The first sum in ([5]) , which we will refer to as the "diago- 
nal part" Pd\t q ), is a statistical mixture of phase states. 
It is mainly responsible for the structure of the phase 
profile given by the Husimi distribution [2l|. The sec- 
ond sum in to be referred below as the "off-diagonal 
part" p^(tq), accounts for quantum correlations and in- 
terference effects, such as, for example, fringes in the 
eigenvalue probability distributions of J x and J y |21[. 



IV. DYNAMICS IN THE PRESENCE OF PHASE 
NOISE 



III. DYNAMICS IN THE ABSENCE OF NOISE 

In the absence of noise, let us describe the quenched 
dynamics of the BJJ induced by a sudden switch off of 
the Josephson coupling K in (pQ) at time t = 0. We take 
A = for simplicity. Initially, the BJJ is in the phase 
state with = (i.e., a = 1). This state is the ground 
state of the Hamiltonian (pQ) in the regime KN ^> x 
where tunelling dominates interactions. In the internal 
BJJ, it can be produced by applying a short 7r/2-pulse to 
the atoms initially in the lower level. Under the effect of 
the quench the phase starts diffusing along the equator 
of the Bloch sphere. The visibility of Ramsey fringes [HI 



We now account for the effect of noise by considering 
the Hamiltonian 

H(t) = X Jz ~ Kt)Jz (6) 

where \(t) is a classical stochastic process. Since 
[H(i),J z ] = at all times, J z is conserved as in the 
noiseless case. Neglecting the fluctuations of Ui (which 
seems justified in the experiments), the fluctuations of 
A are equal to those of AE and are independent of N. 
For a given realization of the process A, the Schrodinger- 
evolved state is obtained from the state \i/j(°\t)) = 
e~ lxJzt \a = 1) in the absence of noise through a rigid 
rotation around the z-axis by a random angle <p(t) = 



-J^dr\(r), i.e., \i/j(t)) = e~^^ |^ (0) (*))- The phase 
4> has a distribution /(</>,£) = j dP [A] 5(<fi(i) — (j)) where 
P [A] is the probability distribution of the process A. Av- 
eraging over all realizations of A leads to the density ma- 
trix p(t) = J dP[X] |^(£))(^(£)|. This is the analog of 
tracing out the bath degrees of freedom in models of sys- 
tems coupled to quantum baths. We obtain 

/oo 
#/(^t)e-^^(°)(t)e^ (7) 
-CO 

where p {0) (t) = |^ (0) (*))W> (0) (*)I is the density matrix in 
the absence of noise. By projecting Eq.([7j) over the Fock 
basis we get 

/CO 
#/(^,i)e- i ^"-"')(n|p(°)(t)|n') 
-CO 

= f\n f -n,t)(n\pW(t)\n') (8) 



where f(m,i) = t)e irn< ^ = e im ^ is the 

Fourier transform of £) with respect to (j) and the 
overline denotes the average over the realizations of 
the noise A according to the probability distribution 
P[A]. To be specific, let us consider a gaussian noise. 
Then f(m,t) = e - a2 ( £ ) m2 /2 e -^m^ w j iere the variance 
a 2 (t) is given in terms of the noise correlation function 
h(r - t') = A(t)A(t') - A 2 = AE(r)AE(r f ) - AE* 
by a 2 (t) = JldrJldT'h(T-r') = 2 J Q £ dr JJ" duh(u). 
(Note that h depends on the time difference r — r' 
by the stationarity of the process, which also implies 
A(£) = A(0) = A.) This yields 



(n\p(t)\n') 



iXt ( n - n '\n\p^\t)\n'). (9) 



The effect of the noise is to spread the noiseless evolu- 
tion (£) along the equator of the Bloch sphere by the 
amount a(£); in this sense, it is a pure-dephasing noise. 
Since our result (|9]) does neither rely on a perturbative 
approach nor on a Markov approximation, it is valid also 
for strong noise and at short times. The variance a 2 (£) 
does not depend on N and completely characterizes the 
effect of the noise on the BJJ. It is given by 

2 ^ \h(0)t 2 if £ < £ c (small time) 

a W ~ [2 J °° drh(r) t if t > T c (Markov) 

where we have introduced the noise time scales t c and T c , 
t c being the largest time such that h(r) ~ h(0) = SX(0) 2 
for |r| < t c and T c the smallest time such that h(r) c± 
for r > T c . Eq. (|9|) shows how the noise suppresses the 
off-diagonal elements of the density matrix in the Fock 
basis. At long times t ^> ( J °° dr/i(r)) _1 the state of the 
BJJ converges to a statistical mixture of Fock states with 
the same probabilities as the initial state, 

N/2 f 27V d6 7T 7T 

n=-N/2 J ° 




OO 
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FIG. 1: (Color online) Phase relaxation of the q = 4 and q = 2 
macroscopic superpositions in the presence of noise sketched 
along the equator = n/2 of the Bloch sphere. Top panels: 
q = 4 (U = T/8) and a 4 = 0, 0.64, 2.05 (from left to right). 
Middle panels: q — 2(ti — T/4) for the same noise intensities 
J °° drh(r) in the Markov regime {a<2 — 0, 0.9, 2.9). The circle 
sizes illustrate qualitatively the phase distribution /(</>, £2,4). 
For intermediate noise (middle column), the superposition is 
closer to the steady state (last column) for q — 4 than for 
q — 2. Bottom panels: Husimi distribution Q(0 = 7r/2, (j)) for 
q — 2 for the same values of a^. Here A = and N = 10. 



The last equality is obtained from ([2|). It means that at 
large times the phase <p is uniformly spread on [0, 27r], as 
is the case for Fock states (FigJTJ right panels). 

For a non-gaussian noise, the two exponentials in the 
right-hand side of (j9]) coincide with the cumulant expan- 
sion of /(m, £) = e im ^^> up to the second cumulant [22] 



and higher cumulants yields extra factors e M^)( n - n ) p ? 
p = 3,4,.... For small times £ < £ c , these factors are 
close to unity and the right-hand side of (|9|) still gives a 
good approximation of the matrix elements of p(t). 

Under the effect of the noise, the visibility (j3]) acquires 
an additional decaying factor due to the above-mentioned 
phase spreading. Indeed, one easily obtains from (j9|) 



(12) 



As found e.g. in superconducting circuits [23| and in 
quantum dots [24], the dephasing factor e _a2 ^/ 2 displays 
a Gaussian decay at short times £ < £ c , corresponding to 
the universal regime of Ref. [HI , and an exponential de- 
cay at long times £ ^> T c (Markov regime). A Gaussian 
decay of the visibility ([T2]) has been observed experimen- 
tally in the internal BJJ at small values of the interac- 
tions x |23| . This indicates that the experiment was per- 
formed in the small-time, non-mar kovian regime. The 
effect of the noise on the visibility decay in this regime 
is shown in FigJ21 for experimentally relevant parame- 
ters 25]. A direct comparison with the experiment re- 



(11) quires to include in the model atom losses. 



FIG. 2: (Color online) Visibility v{t) as a function of time 
(in units of seconds) for x — n ' 0.05Hz, tv • 0.13Hz, tv • 0.25Hz 
(from top to bottom), N = 400. Solid lines: decay of v{t) in 
Ea.(H2l) in the limit \t < 1 and At < 1 with a 2 (t) = h(0)t 2 
and /i(O) 1 ^ 2 = 8Hz. Dashed lines: decay of z/°)(£) under the 
unitary evolution only. For small values of the interactions 
the decay is mainly due to the noise. 



FIG. 3: (Color online) Top: relaxation of the diagonal part of 
the density matrix in the Fock basis (|13[) for q — 2 and N = 10 
to the diagonal matrix as the noise is increased from 

a<i — (left) to a<i — 0.9 (middle) and a<i — 2.9 (right panel). 
Bottom: off-diagonal part of the density matrix for q = 2 
and the same values of a 2 , its vanishing indicates decoherence 
among the components of the macroscopic superposition. 



V. DECOHERENCE AND RELAXATION OF 
MACROSCOPIC SUPERPOSITIONS 

Let us then study the impact of decoherence and relax- 
ation on a macroscopic superposition at the fixed time t q 
as a function of the noise intensity. We assume that q is 
even. Using ((9]) and the decomposition ([5j) of the density 
matrix into diagonal and off-diagonal parts, and setting 
a q = a(t q ), one gets 

<n|p d ,od(t,)|n') = e- a ^^{n\pf oA {t q )\n') (13) 

up to a phase factor irrelevant for decoherence. Since the 
effect of the noise in Eq. ([T3j) factorizes out let us concen- 
trate on the structure of the density matrix in the absence 
of noise. By ((4]) and the definition following Eq.([5j) we 
have {n\pf{t q )\n') oc ££o e^C*-*') = if n' + n 
modulo q and 



(^ 0) (^)K) = ^( n 



TV 



1/2 
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otherwise. Therefore 



a(0 ) , M , x f (n\p(0)\n') if n' = n modulo q 
(n\p\\t q )\n) = < 



if n' 7^ n modulo q 



and 



(n\p%(t q )\n f ) 



(14) 



if n' = n modulo q 

^ (n ' 2 - n2) (n|/3(0)|n') if ri ^ n modulo q 



(15) 

The last equality follows from ([T4|) and from the iden- 
tity p(°\t q ) = e~ ltqXj z p(0)e ltqXj z . In particular, we get 

(n\p^d(tq)\n ± 1) ^ from ([T5]) (since q > 2 is even). 

Thus from ([T3]) - ([T5]) we obtain that Pd(t q ) — >> p(oo) 
when a g g ^> 1 (phase relaxation) and p Q d(t q ) — > when 



a g ^> 1 (decoherence). Hence, in the strong noise limit 
the diagonal part of p(t) relaxes to the steady state (fTTj) 
and the off-diagonal part is washed away (see Figj3j right 
panels). Remarkably, the decoherence factor in Eq. (fT3|) 
does not depend on the atom number N. Also note the 
different noise scales relevant for decoherence, a qi and 
phase relaxation, a q q. When increasing the noise inten- 
sity, Pd(t q ) approaches p(oo) before p d(t q ) vanishes. The 
higher the number of components q in the superposition, 
the more pronounced is this effect. In fact, superposi- 
tions with higher q are better protected against deco- 
herence since they are formed at shorter times and a(t) 
increases with time. Moreover, in the Markov regime 
t q ^> T c , for a fixed noise intensity phase relaxation has 
a stronger effect on states with higher g, as illustrated 
in Fig. CD Indeed, by (HDD, CE}, and (Q~4j), the n ^ n' 
matrix elements of pdifq) are damped by a factor equal 
to or smaller than e~ a< i q I 2 c± exp[— (irq/x) Jq° drh(r)]. 
In the small-time regime t q < £ c , all the ^-component su- 
perpositions relax to p(oo) at the same (g-independant) 
noise intensity (since a q q is independent of q). As a con- 
sequence of the distinct noise scales for decoherence and 
phase relaxation, the BJJ does not turn into a statistical 
mixture of phase states but relaxes directly to the mix- 
ture of Fock states ([TT]) . This is illustrated in Fig J3] for 
the two-component superposition. 

Phase relaxation can be represented by the 
distribution of Pd(t q ) [21(, which for 
TV > 1, and a 2 > iV" 1/4 is given by 

7rA/(2x),e" 2a 2) with 6 3 the Theta 

' 1+sinJ 



Husimi 

q = 2, 
Qd(M) 

Qoo(0)e 3 (- 

function 126 



and Qoo{0) — (- 



the distribution of the state ([TT]) ; Q(7r/2,0) is plotted 
for various values of a 2 in FigJTJ In the absence of 
noise it shows peaks at (j) = and 7r, which correspond 
to the two coherent states of the superposition. The 
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peaks are smeared at increasing <22, and finally at 
CL2 ^> 1 the Husimi distribution reaches the flat profile 

g(7T/2,0) = Q co (7T/2). 

An important consequence of the anomalous decoher- 
ence found in this model is that at intermediate noise 
strength associated to an already important phase relax- 
ation, i.e., such that Pd(t q ) is close to p(oo) (e.g. for <22 = 
0.9 in Figj3]), the system still displays quantum correla- 
tions, which could be exploited in inter ferometry. These 
correlations can be quantified by computing the quantum 
Fisher information Fq[p] associated to the A/"-particle 
density matrix p jHHzj- The quantum Fisher information 
allows to estimate the best possible phase sensitivity AO 
in an interferometric scheme according to the generalized 
uncertainty principle AO > 1 / (Vm^/ Fq\p\) , m being the 
number of measurements performed [27j. For separable 
states we have F Q [p] < N and AO > A0 S n = 1/VmN. 
Sub-shot noise sensitivities AO < AOsn can be achieved 
for states p satisfying Fq[p] > N, which is also a suffi- 
cient condition for multi-particle entanglement [H|. For 
d2 = 0, the two-component superposition p^fo) has 
Fq = TV 2 [111, corresponding to the Heisenberg limit 
AO = l/\/mN. For noise strength <22 = 0.9 and TV = 10 
we find [28, 2§| FQ[p(t2)] — 58, which is indeed still larger 
than FQ[p^\t2)) = N and leads to a sensitivity gain 
AO/ AOs n of —3.8 db with respect to the use of sepa- 
rable states. In the limit a2 ^> 1, pfo) — p(oo) has a 
Fisher information Fq = N(N - 1)/(2N + 2) < N @. 

VI. SUMMARY AND CONCLUDING REMARKS 

We have treated analytically the time evolution of a 
phase state driven simultaneously by the atomic inter- 
actions and by a noise coupled to the number imbalance 
operator J z . We have derived an exact expression for the 
density matrix, as well as for the noise-induced decay of 
the Ramsey visibility. The effect of the noise on the cre- 
ation of macroscopic superpositions of phase states is to 
cause decoherence, i.e., the vanishing of quantum corre- 
lations. We have found that if the fluctuations of the 
atomic interactions are negligible, decoherence is "less 
efficient" than phase relaxation, especially for superpo- 
sitions with a large number of components. As a conse- 
quence, the states generated by the noisy dynamics could 
in principle lead to sub-shot noise precision in interfer- 
ometry. 

The surprising fact that decoherence is not enhanced 
by increasing the atom number N is specific to the noise 
considered. Indeed, such a noise is applied perpendicu- 
larly to the equator of the Bloch sphere where the phase 
states of the superpositions lay. As a result, the noise is 
insensitive to the separation between these states, which 
scales with N. However, such superpositions are very 
fragile under a noise applied parallel to the equatorial 
plane, which resolves the separation between the compo- 



nents. This yields an indication as to which classical noise 
to reduce to preserve the coherence in superpositions of 
the phase states: this is the noise in directions parallel to 
this plane. For example, stochastic fluctuations on the 
tunnel amplitude K give rise to rapid decoherence of the 
macroscopic superposition (\a = 1) + e %1 \ot = — l))/y/2 
at a rate increasing with the atom number, without in- 
ducing relaxation. 
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